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Abstract 


Visco-eiastic polymers like natural and synthetic rubbers ( with various additives ) are 
commonly used as vibration isolators. Such rubber elements provide both stiffness and 
damping as required by a good isolator. However, proper design can only be made only if 
their dynamic characteristics are correctly modelled. 

In this thesis, a hollow cylindrical rubber specimen is tested under harmonic excitation 
for modelling its dynamic characteristics. A dynamic constitutive law, under harmonic 
loading, is proposed to fit the experimental results. Both the amplitude of deformation and 
frequency are varied in the range of interest for mechanical structural applications. Marked 
non-linearity is exhibited. The dependence of dynamic stiffness on the amplitude and a 
combination of linear hysteretic plus nonlinear Coulomb damping are seen to be essential 
for good modelling. For satisfactory modelling, future work should al^o take into account 
the effects of temperature. 
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Chapter 1 


Introduction 

Resonance induced vibrations are detrimental in most of the engineering situations and 
may even lead to drastic failures. The sustenance of vibrations over long periods may 
cause sudden fatigue failure. The need for control of vibrations also arises in high precision 
measuring equipments. Frequently, one encounters the problem of protecting a machine 
or equipment from the vibrating foundation. On the other hand, the foundation may 
sometimes need to be protected from the vibrating machinery. An effective way of curbing 
the transfer of vibrations from the vibrating source is to have a resilient element with or 
without a damping element, known as an isolator. This method of vibration control is often 
used if it is not possible to eliminate the vibrations of the source. 

Sometimes an isolator is designed to have both the resilience and damping. With re- 
markable upswing in the rubber technology, various synthetic rubbers .and polymers ha.ve 
come up as promising materials for isolators as it is possible to combine both the proper- 
ties in a single piece, unlike the structural materials where it is necessary to have separate 
dampers or inserts of high energy dissipation capacity materials since the damping capacity 
of structural materials is usually very small. 

Rubber offers several other useful properties making it suitable to a large class of engi- 
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neering problems. It has extensive use in fatigue and noise control in industry e.g. drilling 
rods, grinders, circular saw blades. Light weight combined with above-mentioned qualities 
enhances its adoptability in the applications having weight reduction as a \'ita] and un- 
comprisable objective such as aircraft industries or space applicatins. Its combination with 
structural materials is also in vogue. 

Rubber and polymer properties undergo considerable change by variation of their con- 
stituents and hence desirable properties can be achieved. These properties are heavily de- 
pendent on temperature and frequency. To accommodate both temperature and frequency 
effects in a rubber or polymeric material, temperature effects are reduced to equivalent 
frequency with the help of master curves [1, 2] but even at a constant temperature sev- 
eral factors affecting the properties of these materials, namely, frequency, amplitude, time 
history, aging and type of excitation, demand attention from the designer. Application 
of rubber and polymers in a particular problem needs indentification of requirements of 
the problem and then judicious selection of material, shape, geometry and dimensions of 
isolators, absorbers or any other vibration controlling element which needs prediction of 
the behaviour of these controlling elements in different sets of operating conditions [3, 4]. 
consequently, the need for the precise models for these materials accommodating aforesaid 
effects becomes manyfold. 

Nevertheless, there has been greater motivation toward the study of rubber and poly- 
meric materials in recent years which is reflected in the advances in the phenomenological 
theory of linear visco-elasticity and spectacular success in the synthesis of rubbers for at- 
tainment of stringent requirements imposed by specific problems. 

The approach to the design changes for different applications. The modelling is con- 
cerned with obtaining a mathematical statement of the properties of the system, best ex-- 
pressed through a rheological constitutive equation, correlating the input and the response. 



To avoid undue complications in the model, some simplifying assumptions are made and 
the system is idealized. Fixing the domain of the study is important from the view-point 
ol pursuing solely the properties of interest. The analysis of the gross behaviour of the sys- 
tem is carried out by involving the basic unit properties of the component materials. The 
understanding of solid state or micromechanistic concepts for an engineer is also gaining 
significance, despite the fact that this approach has still not flourished to the point where it 
could become the sole basis of the complete design process. Ad-hoc testing can sometimes 
be taken up after analyzing the viability of the problem in terms of time, size and cost. 

While specifically dealing with the modelling of rubber-like materials, its elastic proper- 
ties should be studied with emphasis on the buckling behaviour. Hollow rubber specimens 
are vulnerable to buckling in compressive loading and generally have a low tensile strength. 
Abrasion resistance and other similar properties should also be given due attention if the 
application so requires. The damping in general is the most difficult to model. The tem- 
perature dependence of the elastic or loss modulus should be established, especially, if the 
transition region is near the room temperature. 

1.1 Review of Previous Work 

Recently, Tinker and Kutchins (1992) [5] investigated the damping phenomenon in wire rope 
isolators. They have reported that damping in typical wire rope isolator is predominately 
due to Coulomb friction. They developed a serai-empirical model involving power 
damping and variable Coulomb friction damping. Harris and Stevenson (1986) [6] published 
the results for natural rubber (NR) and NR blends with high damping polymers. NR 
was found to be nearly linear in simple shear in both dynamic as well as static tests. 
Its behaviour was close to that of a steel spring but peak of the transimissiblity curve 
was low due to relatively high damping. Filled NR was drastically iionlinear and showed 



Introduction 


(I(’|)(’IkI(’Iic.o oti frcHjiu'iicy and amplitudn of strain. As an cxtr’iision to this work, Harris 
(1987) [7] investigated the effect of non-sinusidal input on rubber and noted that for complex 
vibrations (■f)nHiHting of tnany amplitudes and freqiieticies, the overall dyiiainic. behaviour 
tends toward that obtained in single sine tests at a strain amplitude which corresponds to 
the largest amplitude that occurs in time history. Under a. more and more complex system 
of vibration, nonlinear rubber shows greater linearity in behaviour and exhibits higher 
damping than indicated by their dynamic properties measured in conventional sinusoidal 
tests. Therefore, when the service enviornment of many components is taken into account, 
the dynamic behaviour of nonlinear rubber may give better performance characteristics 
than expected from conventional sine tests. 

Medalia (1978) [8] discussed the effects of Carbon Black on dynamic properties of rubber 
vulcazinates and proposed empirical relations for loss modulous of rubber. In reference [9], 
elastic modulous and loss modulous of SBR (Styrene Buta-diene Rubber ) containing parti- 
cals of cross-linked polysterene have been obtained. Vinogradov and Pivovarov (1986) [10] 
introduced a hysteretic model comprising a linear viscous part and a nonlinear Davidenkov 
part in order to simulate rate and amplitude dependent internel losses. For a vibrating 
cable, parameter identification has been discussed and theoritical and experimental results 
have been compared. lyenger (1979) [11] suggested an analytical model for the hysteretic 
behaviour of beams. The model was used to simulate several types of realistic softening 
hysteretic loops. 

1.2 Present Work 

This work is aimed at modelling a cylindrical rubber sleeve, commonly known as rubber 
spring, under harmonic loading. The rubber spring was used to isolate a known mass from 
a harmonically vibrating base plate. The objective was to obtain a constitutive equation 
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r, 

correlating the force with the instantaneous deformation and the instantaneous rate of 
deformation of the rubber isolator. For this, the base plate was sinusoidally excited using 
an electrodynainic shaker. The base displacement, top displacement and top acceleration 
signals were recorded and digitized. The deformation i.c. the relative displacement across 
the isolator and hysteresis loops were obtained by using the digitized displacement and 
acceleration data. The variation of area and shape of the hysteresis looj)s with frequency 
and amplitude was studied on the basis of which, a form of the damping was ascertained. 

The frequency range under consideration was 20 Hz to 200 Hz which is also the range of 
interest in most of the industrial applications. The amplitude dependence of rubber proper- 
ties has been established. The behaviour of the proposed model aggrees with experimental 
findings. 

The thesis has been organised in the following way. The first section of Chapter 2 
briefly discusses the general theory behind the modelling of an isolator and commonly 
prevelent methods thereof. The second part of the chapter defines the present problem 
and outUnes the intended method of modelling. In Chapter 3, the description of the set- 
up and the instrumentation is given. The processing of the displacement and acceleration 
signals obtained through the experiments has been explained. Chapter 4 discusses the 
observations of the experiments and inferences drawn on the basis of them. The last section 
of this chapter analyses the inferences and establishes a mathematical model for the isolator. 
Chapter 5 concludes the thesis with a brief discussion about the scope of future work. 



Chapter 2 


Modelling of an Isolator through 

Harmonic Loading 

A system undergoing mechanical vibrations has an input or forcing function and the re- 
sponse as its effect, as illustrated in F’ig. 2.1. The relation between the forcing function and 
response depends on the properties of the system. 



Figure 2.1; Relationship between forcing fnnctiion and response 

The modelling of the system defines the relationship between the input and the output 
in the form of a law or a constitutive equation. Thus, the modelling deals with the deter- 
mination and the mathematical specification of the properties of the system from the input 
and response characteristics. 
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Modelling of an Isolator through Harmonic Loading 7 

2.1 Modelling of an Isolator 

The simplest type of a vibrating system is a single degree of freedom system in v^hich the 
state of the system is completly determined by only one dependent variable. A single- 
degree-of-freedom system is modelled as a resilient member with or without energy dissi- 
pation mechanism connecting a mass and the foundation. An externel force acts on this 
system which sustains the vibrations. The resilient member is known to be an isolator 
when it is used to reduce the transfer of ill-effects of vibrations from the base to the mass 
( an equipment or machinery in practicl situations ) or vice-versa. Two types of problems 
normally encountered in practice, for a harmonically excited single degree of freedom sys- 
tem, are shown in Fig, 2.2(a) and 2.2(b). In one class of the problems, the base on which 


Fo Cos wt 



Moving Base Fixed Base 


(a) (b) 


Figure 2.2: An isolator under vibrations (a) Base Excitation Case (b) External force exci- 
tation case 


the system is mounted undergoes vibrations and causes the system to vibrate while in the 





other, an exte.rne] force causes the vibrations of the system. Both of these problems in a 
linear system are mathematically equivalent. 

A vibrating system may be subjected to several types of loading such as sinusoidal, 
saw-tooth wave force, a general periodic force or a random force. Out of these, sinusoidal 
loading is commonly encountered in the engineering practice as this is generally caused due 
to unbalance in rotating machineries which are widely used in the industry. 

There are three properties of the system, viz. inertia, stiffness, and damping which 
govern the constitutive equation of a vibrating system. The general form of constitutive 
equation for a harmonically excited single degree of freedom system is written as 

rnx F(x, x) = Focoscvt (2-1) 

In the base excitation case [ Fig. 2(a) ], it can be alternatively written as 

mx -f F{z, z) = 0 (2.2) 

where z — x — y is the relative displacement of the mass with respect to the base and 
z = X - yis the relative velocity. The Fiz,z) term involves 

(i) Restoring force Fa{z) which depends on the stiffness of the system. This is a function 

of the deformation of the isolator, z and does not involve z. 

(ii) Damping force term Fd{z,z) which is responsible for the energy dissipation of the 

system. This is rate dependent and may involve both z and z. Iherefore, i^( 2 , i) is 
written as 

F{z,z)=Fd{z,z)^-F,{z) 

Equations 2.1 and 2.2 may, in general, take the form of a nonlinear differential equa- 
tion. This nonlinearity depends on a number of factors. In particular, material of the 
isolator plays an important role. For example, in rubbers and polymers the stiffness is 
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found to he nonlinear.! niaterial.s do not follow Hooke’s law and may show hard- 
ening or softening characteristics with increasing strain. Structural materials have 
linear stiffness in the elastic range but may considerably deviate from linearity in the 
plastic range. The dynamic properties of rubbers and polymers have been observed to 
be heavily dependent on the frequency. The stiffness of these materials may undergo 
large change depending on the frequency. It may or may not depend on the stress 
amplitude. Similarly, the damping which is a measure of energy dissipation capacity, 
may also depend on the frequency and amplitude of vibration. If either damping or 
stiffness depends on the vibration amplitude, the material shows nonlinear behaviour. 

Based on the above-mentioned characteristics, damping may be of three types. 

(a) Rate independent or frequency independent which is also called the hysteretic 
damping. 

(b) Linear dashpot damping or viscous damping in which the damping force increases 
linearly with frequency. 

(c) General frequency dependent damping which is observed in many anelastic and 
visco-elastic materials. In this, the damping energy is a function of frequency. 
Case(b) is only a special case of this. 

In linear damping, the energy dissipation is ])roportional to the square of the am- 
plitude. Various synthetic rubbers, polymers and rubber-like materials exhibit the 
third type of damping. For these materials, the force-deformation relationship is 
rate-dependent. They arc modelled by a combination of viscous and elastic elements 
in order to represent both datnping and stiffness. In the most general case, the force— ■ 



deformation relationship in a visco-elastic material is expressed as 


ao + ffli 


\dtj 


+ 0,2 


d^' 

dt^ 


+ 


a = 


bo 4 - bi 




(2.4) 


For harmonic excitation a — croe^^^ and c = this reduces the above equation 

to 

l^ao -f oi (joj) -f 02 + ■■•] Oo = [ 6 o + bj (Jcv) -f 62 (Joj f + ...] €0 (2.5) 


or 

E* = — =z [^l(‘^) + 

Co [Ax{uj) + jA2{ijj)] 

or 

E* = E' {02) ^jE"{u) 

This can be alternatively expressed as F = ICx or 


( 2 . 6 ) 


(2.7) 


F ={K'{o:) + jK"{u:))x ( 2 . 8 ) 

K* is known as the complex stiifness, K'{oj) is known as the storage modnhis and 
associates force F in phase with the displacement x. K"(u) is known as the loss 
modulus and associates the force 90° out of phase with the displacement x. The 
visco-elastic properties show remarkable dependence on the temperature also, so K' 
and A'" may be functions of temperature as well. 

The major shortcoming of the above model with a’s and 6’s as constants is that 
they cannot accomodate stress and strain dependent behaviour observed in several 
nonlinear materials. The properties of some of these nonlinear materials are affected 
by the stress-history, temperature, static mean stress, frequency or the amplitude of 
the stress. 

If F(z, z) is plotted against z [ refer equation 2.2 ], the curve obtained is a loop 
which is called hysteresis loop. The shape and area enclosed by the loop can reveal 


nuich about the above-mentioned characteristics and consequently, are of paramount 
importance in the modelling of nonlinear materials. The area of the loop gives the 
energy dissipated per cycle. The damping energy per cycle can be expressed as 

^ Fd{z,z)dz (2.9) 

The linear damping exhibits elliptical hysteresis loops whose area is proportional to the 
square of the deformation amplitude. In other words the damping energy is expressed 
Eiamp 0^ Z”' with n = 2 and Z is the amplitude. In some polymers and elastomers, 
n is greater than 2, but the shape of the loop remains close to elliptical. In others, 
not only is n greater than 2, but the shape of the loop is very different from elliptical. 

One imporatant class of nonlinear damping is Coulomb Friction damping in which 
the damping force is independent of the magnitude of the velocity. For this type of 
damping, the energy disssipation is proportional to the amplitude i.e. n = 1. The 
hysteresis loop obtained is a parallelogram if the stiffness is linear. 

The inclination of the hysteresis loop indicates the stiffness of the system and is 
governed by Fs{z) [ refer eqn. 2.3 ] whereas the area of the loop is governed by 
Fi{z,z). Figures 2.3, 2.4 and 2.5 explain how the hysteresis loop in Fig. 2.5 is 
obtained by superimposing Fa(z) vft. z curve in Fig. 2.3 and Fd(z,z) tw. z, shown in 
Fig. 2.4. 

Ilcsidos energy dissipation per cycle i.c. the area of the hysteresis loop, other measures ., 
of damping capacity are sometimes used. Some of them are less common and are 
applicable to special class of problems only. Some of these measures are logrithmic 
decreement, half-power band width and resonance magnification factor. The loss 
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Force 



Figure 2.3: Elastic restoring force F* vs. z 


factor is a widely used measure of damping. It is defined as 

D 

~ 2ttU 



( 2 . 10 ) 


where D is the damping energy per cycle and U is the maximum strain energy stored 
in the system during the cycle. The dependence of damping energy, stiffness and other 
properties on frequency, temperature, amplitude, mean static stress etc. is required 
to be considered for the modelling problem. 


2.2 Present Problem 


The objective is to model a rubber isolator. The isolator was made of natural rubber 
with Carbon Black fillers. The rubber isolater was mounted on an electrodynamic 
shaker and was fixed with a mass at the top. The base, on which the isolator was 
mounted, was harmonicafiy excited. The sinusoidal excitation was selected for the 
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Force 



Figure 2.4: Dissipative force Fd vs. z 

experiments because sinusoidal loading is encountered in a number of engeermg prob- 
lems. Secondly, non-sinusoidal periodic loading can be represented as the summation 
of its Fourier components and thus can be analyzed through superposition of the ef- 
fects of individual sine components. Therefore, this system can be modelled as a single 
degree of freedom system harmonically excited at the base. The freebody diagram 
of the system is shown in Fig. 2.6. The rubber isolator is attached to the moving 
base and the force acting on the rubber isolator is expressed as a function of the rel- 
ative displacement and the relative velocity of the mass with respect to the base i.e. 
as a function of the deformation and rate of deformation of the rubber spring. The 

constitutive equation from equation 2.2 can be written as 


mx + F{z,z) = 0 
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Figure 2.5: Combined effect F{z^z) vs. -ar 

In order to determine the exact form of F(< 2 r, i), acceleration x of the mass and defor- 
mation z across the isolator are needed. Two pick-ups,- one at the base and the other 
at the top, were used for obtaining x and y. The hysteresis loop is obtained by 
plotting F(z^z) against where F{z^z) is given by —rnx. The area of the hysteresis 
loop was used as a measure of the damping capacity, experiments were conducted for 
studying (i) the damping energy with frequency behaviour keeping base displacements 
such that the relative displacement across the isolator i.e, the deformation, is constant 
(ii) the damping energy with deformation of the isolator at constant frequencies. 

Although rubbers and polymers may have large dependence on the temperature, the 
temperature dependence was not undertaken for the study and the experimentas were 
conducted at room temperature. Similarly, the stress-history effect was also not in- 
vestigated. 
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Figure 2.6: Freebody diagram of the isolator 





Chapter 3 


Experimentation 


3.1 Set-up 

A cylindrical sleeve made of Natural Rubber with Carbon Black fillers ( IRHD 40 ) 
was used as the isolator. The external and internal diameters of the rubber spring 
were 34.8 mm and 21.2 mm respectively and the height was 98.1 mm. The sketch of 
the rig is shown in Fig. 3.1. The rubber cylinder at one end was attached to a cast 
iron plate (2) through hose clamp and adhesive. The rubber sleeve-cast iron plate 
assembly was fixed in a housing (3) serving as the support for the rubber sleeve and 
the top plate. The projected ends of the base of the housing have two open U-slots. 
A base plate (5) was constructed for facilitating the mounting of this assembly on an 
electrodynaraic shaker. The base plate has seven countersunk through-holes for fixing 
it at the to|) of the shaker. 1’he housing cf)ntainiiig the rubber sleeve was mounted on 
the base plate by fastening two bolts ( not shown in the figure ) through U-slots of 
the housing and two additional holes in the base plate. In order to increase the mass 
at the top of the rubber cylinder, a mild steel block (1) was screwed to the top of the 
plate (2). 
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3.2 Instrumentation 


As (liscussnd in (Jhapter 2, acceleration x of the mass and the deformation z of the 
rubber spring are needed for obtaining the hysteresis loops. For this, three signals 
viz. acceleration and displacement at the top and the displacement at the base, were 
recorded. Two B & K 4370 accelerometers, one at the top plate and ther other at the 
base plate, were used to pick up the acceleration signals as shown in Fig. 3.2. These 
signals were amplified by using two B & K 2635 charge amplifiers and integrated 
twice to give the displacement signals. The deformation of the isolator was obtained 
by taking the difference of the displacement signals at the top and the base. Since 
the phase of the signal may undergo some shift in the amplification and integration 


[)rocess, it was ensured that both the amplifiers were always at the same settings which 
would reduce the possibility of unequal phase phase shift in the two signals, because 


any error in the relative phase of the two signals could lead to erroneous results. The 
lower and higher cut-off frequenies were 1 Hz and 3 KHz in the displacement mode and 


2 Hz and 3 KHz in the acceleration mode. The. amplification in the displacement mode 
was 1000 mV/unit out till 120 Hz and 100 mV/unit out at higher frequencies. The 
amplification in the acceleration mode was 10 mV/unit out for the entire frequency 

range. 


The signals were observed on Kikusui COR 5501 U Two Channel Digital Storage Os- 
cilloscope. The oscilloscope had a 4 kilobyte memory for each channel, so the signals 
could be stored simultaneously. The signals were digitized and transferred to a mi- 
crocomputer through RS232C. A waveform on the CRF screen c)f tlio. oscilloscope bad 
4096 (horizontal) x 256 (vertical) resolution. Therefore, a digitized signal could have 
4096 data points out of which 4000 points were recorded. 
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Figure 3.2: Instrumentation set-up 

The acceleration of the top plate was also to be recorded along with the displace- 
ments at the top and the base. However, since the oscilloscope had only 2 channels, 
all the signals could not be recorded simultaneously. The two displacements signals 
were recorded separately. But this caused the incompatibility in the phase of the ac- 
celeration signal and the two displacement signals as they were recorded at different 
time instants. The incompatibility of the phase was corrected later during the signal 
processing which is discussed in the next section. 

The base plate was excited by I. M. V. electrodynamic shaker having pneumatic 
suspension base. The frequency was monitored by the digital display on the control 
panel of the shaker and cross-checked with the oscilloscope using vertical cursors on 
the CRT screen of the oscilloscope. 

The damping energy variations (i) with frequency ( from 20 Hz to 200 Hz in steps 
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of 20 Hz ) at a constant deforniation amplitude of the rubber spring, and (ii) with 
deformation amplitude ( 0.05 mm to 0.3 mm at every 0.05 mm ) at a constant fre- 
quency, were studied. In the first set of experiments, three values of the deformation 
amplitude, namely, 0.05 mm, 0.1 mm and 0.15 mm, were used to study the damping 
energy with frequency behaviour. In the latter catagory, three representative frequen- 
cies, viz. 60 Hz, 80 Hz and 120 Hz were used. The amplitude of the deformation was 
monitored on the oscilloscope using subtraction mode of the oscilloscope which gives 
the difference of the signals being supplied to the two channels. 

A list of specifications of the equipments used is given in the appendix-A. 

3.3 Signal Processing 

After digitizing the three signals, the high frequency noise was filtered out using linear 
phase-shift, causing equal phase shift in all the signals, low pass filter. The cut-off 
frequency used was 3 kHz. A signal before and after filtering is shown in Fig. 3.3. 



It has been discussed in Chapter 2 that the hysteresis loop is obtained by plotting 
p~'{z,z) against deformation z. F(z,z) equals —mx which means that it is 180° out 
of |)hasc with the acceleration x of the mass at the top. Acceleration x in turn is 
180” out of phase with the displacement x of the top, thus the force F{z,z) and the 
displacement x are in the same phase. Therefore, the recorded acceleration signal was 
processed to have the same phase as that of the displacement signal of the top plate 
and multiplied by mass m, so as to give a measure of the force F{z,z). Since the 
acceleration signal at the top was not measured simultaneously with the displacement 
signals, the initial phase of the acceleration signal ( Fig. 3.4 ) is synchronized with 
that of the displacement signal. 

After this, the deformation of the rubber was obtained by subtracting the base dis- 
placement from the top displacement. F{z, z) obtained from the recorded acceleration 
by the above procedure was plotted against z to get the hysteresis loop. Its area was 
calculated numerically using trapezoidal rule on the data points recorded in the PC. 
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Figure 3.3: An example of filtering ( base displacement signal of 0.15 mm amplitude at 120 





Chapter 4 


Results and Discussion 


4.1 Static Stiffness 

As discussed in chapter 2, the modelling of the system requires the determination of 
the term F{z,z} which consists of the damping term Fd(z,z) and the stiffness term 
Ft{2), The stiffness is one of the three properties of the system goveraing its behaviour 
in vibrations. In a large class of rubber-like materials, the stiffness of the material 
in dynamic testing may show marked departure from the static stiffness, yet, static 
stiffness gives an overall idea of the nature of the force-deformation relationship of 
the material. 

The static stiffness test in compression was conducted on Universal Testing Machine 
INSTEON 1195. The maximum comprtession was limited to 0.5 mm since in the dy- 
namic test, the maximun deformation amplitude was 0.3 mm. It should be emphasized 
that the static stiffness experiment could not be done accurately because the sidewise 
deformation of the rubber spring could not be entirely eliminated. This resulted in 
some eccentricity of the loa<iing, which violated the condition of nniatriai stress in the 
static stiffness test in compression. The experimental static stiffness curve is shown 
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Figure 4.1: Static stiffcess curve 

in Fig. 4.1. It is observed that the rubber spring is softening type. 

4.2 Natural Prequency 

The natural frequency of this system can be qnite accurately determined by deter- 
mining the frequency at which the relative transmissibility of the system reaches the 
maximnin valne for a constant value of the base displacement amplitude. The spring- 
mass system was excited with a constant base amplitude of 0.1 mm and the relative 
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Figure 4.2: Natural frequency curve 

amplitude of the mass i.e. the deformation amplitude across the rubber spring was 
measured at different frequencies. The deformation amplitude is plotted against the 
frequency in Fig. 4.2. It is observed that this response curve is typical of the Voigt 
model. The natural frequency from Fig. 4.2 is observed to be 57 Hz. 




4.3 Frequency Dependence 


The. frnqiio.nry (Inpoiulonro of thn HUfFnofiR of tlio nihbor spring rnal.orial is host (rb- 
served through the hysteresis loops. Fig. 4.3 shows three hysteresis loops obtained 
through the procedure described in section 2.2 and 3.3. The loops 1, 2 and 3 are at 
constant frequencies 50 Hz, 90 Hz and 150 Hz with the base excitation being set up 
to such amplitude that the maximum relative displacement across the rubber spring 
i.e. the deformation of the rubber spring in each case is obtained as 0.1 mm. These 
hysteresis loops resemble ellipses but are pointed at the ends. Largely elliptical shape 
of the loops indicates the dominance of the viscous or power damping. The in- 
clination of the loop which is governed by the stiflhess of the system remains nearly 
unchanged with frequency. The constant inclination of the hysteresis loops was ob- 
served at other values of the deformation amplitudes also. 

The area of the hysteresis loops gives the damping energy of the specimen per cycle. 
The damping energy variation with frequency at constant deformation amplitude 
are plotted in figures 4.4 and 4.5 for 0.05 mm and 0.1 mm deformation amplitude, 
respectively. It is seen 

that there is no clear trend in the behaviour of the damping energy variation with 
frequency. However, the scatter in Figs. 4.4 and 4.5 can be seen to be about a mean 
value of constant energy dissi|)ation. Therefore, in the absence of a clearctit trend in 
Figs. 4.4 and 4.5, it is reasonable to assume that the energy dissipation is independent 
of frequency in the small range ol the frequency considered in this work. 
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Figure 4.3: Hysteresis loops ( 0.1 mm deformation amplitude; 50 Hz, 90 Hz and 150 Hz) 
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frequency ( Hz ) 


Figure 4.4: Damping Energy vs. Frequency ( 0.05 mm deformation amplitude ) 



Figure 4.5: Damping Energy vs. Frequency ( 0.1 mm deformation amplitude ) 
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4.4 Amplitude Dependence 

In order to study the amplitude dependence of the system, experiments were con- 
ducted for six different values of deformation amplitude of the rubber spring ranging 
from 0.05 mm to 0.3 mm in steps of 0.05 mm, keeping the frequency constant. The 
energy dissipated per cycle is plotted against the deformation amplitude for 60 Hz 
and 80 Hz in figures 4.6 and 4.7, respectively. 

The solid line in the figure is best fit line. The slope of this line is p assuming 
Edamp « \Eo\^ where \Zo\ is the deformation amplitude. It can be seen that at both 
the frequencies, the damping energy variation with amplitude has an index, p, less than 
2 which can be calculated by The departure of p from 2 implies that this 

system cannot have the full contribution of linear viscous or hysteretic damping which 
have an index of 2, as discussed in Chapter 2. In power damping, the damping 
force is given by Fd = (|2:i"spn(3;)). This gives the damping energy Edamp oc 
Thus with p less than 2 for the rubber spring, as observed above, the index n should 
be less than 1. 

Figure 4.8 shows the hysteresis loops at 60 Hz for five different values of the deforma- 
tion amplitude. An important observation from Fig. 4.8 is that the inclination of the 
loop decreases with increasing deformation amplitude. The dynamic stiffness curve 
[12,13] or the backbone curve is the curve obtained by joining the midpoints of those 
points on the upper and lower branches of the hyseresis loop that have same values 
of the displacement. 

Figure 4.8 exhibits that all the hysteresis loops closely resemble ellipses. This indicates 
that the backbone curves are straight lines along the major axes. Similar tests were 
conducted at 80 Hz and 120 Hz and again the shape of the hysteresis loops was 
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Figure 4.7: Damping energy vs. Relative amplitude ( 80 Hz, Solid line is the best fit line) 
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nearly elliplicaJ. 'J’lie fact tliat llie shape of llie hysteresis l()o|) does not undergo any 
noticeable change also indicates that the form of the damping force remains same with 
changing amplitude and frecpiency. 

4.5 Proposed Model 

The main objective of the present work is to frame the constitutive equation for the 
rubber spring which gives an idea about the stiffness properties of the rubber isolator 
and exact damping mechanisms responsible for its dynamic behaviour. 

In order to formulate a model for the rubber spring, the observations from figures 4.1 
to 4.8 ( discussed in earlier sections ) are reproduced briefly here. 

Ob. 1 Static stiffness test ( Fig. 4.1 ) indicates a softening type spring. 

Ob. 2 Dynamic stiffness is indipendent of the frequency ( Fig. 4.3 ). 

Ob. 3 Dynamic stiffness decrea.ses with increasing deformation amplitude Zq ( Fig. 
4.8 ). 

Ob. 4 The dynamic stiffness is linear with the deformation z at a particular defor- 
mation amplitude Zo ( Figures 4.3 and 4.8 ). 

Ob. 5 The damping energy index is less than 2 ( Figures 4.6 and 4.7 ). 

Ob. 6 The damping energy is independent of the frequency ( Figures 4.4 and 4.5 ) 
in the range of frequency normally encountered in the structural vibrations. 

Ob. 7 The contribution of different damping mechanisms is independent of frequency 
and the deformation amplitude ( Fig. 4.8 ). 

The observations (1) and (3) rule out the possibility of using stiffness from the 
static test because F,{z\ modelled with the sialic stiifness, cannot take care of 
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Figure 4.9: Components of F{z,z) 

the rotation of the loops with increasing amplitudes. It is worthwhile to recall 
here that the model F{z,z) ( eqn. 2.3 ) consists of Fs(z) and Frf(2,i), as is 
illustrated in Fig. 4.9. It is apparent that only the term Fg(z) can account for 
observation (3). Now observation (4) suggests a model F,(z) = KdynZ where 
f'^dyn should show dependence upon deformation amplitude Zo ( Observation (3) 
). Therefore, lidyn was selected as /(jj,,! — a — where a, b and g will be 

determined through least square fit. 

The damping model Fd{z,z) is taken to be consisting of two parts, one part rep- 
resenting hysteretic damping and the other part representing Coulomb damping, 
because of the following reasons. The elliptical shape of the hysteresis loops along 
with the observation (6) directly indicates the presence of hysteretic damping. 
However, the observation (5) indicates that in conjunction with the hysteretic 
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damping ( for which p = 2 ), there should be another damping term having 
the specific damping energy exponent p less than 2. Since Coulomb damping 
yields the damping energy exponent p = 1, the damping model is taken to be a 
combination of hysteretic and Coulomb damping such that 

Fd{z,z) = fcoisgn{z)-\- z (4.1) 


However, the above model needs another rerinement because eqn 4.1 will result in 
a hysteresis loop with sharp corners in contrast to the hysteresis loops observed 
in Fig. 4.8. The sharp corners result because of jump in fed at i = 0. 

This is noteworthy that the dry friction fcoi being independent of the 

magnitude of the velocity is the idealization of the actual Coulomb friction be- 
haviour [14, 15]. Due to different contributions of the molecular friction and 
effects of the fillers in the rubber, the form of the Coulomb friction may vary 
with frequency and amplitude, and this is distinctly dilfcrent from the structural 
materials [16]. In the proposed model, too, the Coulomb friction is modelled as 
increasing with velocity wihtin a small velocity zone starting from zero at zero 
velocity and dropping to a constant value with further increase in the velocity. 
This kind of behaviour of the dry friction has been reported in reference [17] 
where the dry friction characteristics of polymeric seals have been discussed. 
Therefore, finally expressed s\h 


Fd = (^) [/co< + - [froi + sgr^z) (4.2) 

The behaviour of the. Coulomb frir.Uoii in tlii.s model is shown in Fig. 4.10. In 
order to confirm to the observations (5) and (6) of the damping energy beliaviour, 
and h in the above term should not be much different from those in equation 
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4.1. Also Cl, C 2 , C 3 and C4 are so selected ( through trial and error ) that 
these constants do not affect the energy dissipation over a cycle but at the same 
time, do contribute in rounding the corners in hysteresis loops. With the above 
observation, tfie energy dissipation per cycle can be written as 

^damp ~ ^h|^o| d" ^fc.ol\^o\ (d*d) 

The constitutive equation of the system is expressed as 

^+[fcol + ^ sgn(z)+K^yr,z = 

(4.4) 

where 

I<dyn = fl - ( 4 - 5 ) 

m, the mass of the plate and the steel block on the rubber spring including the 
mass of the accelerometer, was measured to be 2.2 Kg. For constants a, 6 , q 
in equation 4.5, the dynamic stiffness Kdyn at 60 Hz for different deformation 
amplitudes \Zo\ was calculated ( Fig. 4.8 ). Thereafter, the method of least 
squares was used for finding out o, b, q in equation 4.5. 

From equation 4.3, the damping energy Edamp of th® model ( equation 4.4 ) is 
given as 7 r/i|Zop + 4fcoi\Zo\. The constants h and fcol were determined from 
the damping energy variation with deformation amplitude \Zo\ at 60 Hz ( Fig. 
4.6 ) by least square fit. The damping energy index in Fig. 4.6 is 1.5 and that 
obtained by calculated values of h and fcoi is 1.48. 

Constants C'l, C 2 , C 3 , C 4 were calculated by trial and error. It is observed that 
C 3 gives the value of the velocity around which the Coulomb friction term attains 
the maximum value and this should be small so as to lie near zero. T hough both 
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C 2 C 4 influence the variation of friction with velocity, yet, C 2 dominates the 
behaviour of the friction force after attaining the maximum value and Ca affects 
its increase from zero to the maximum value. The values of the above constants 
are tabulated in the table 4.1. 

Since this model involves wliich is a function of the deformation amplitude 
\Zo\^ it has to be solved by trial and error. The experimental and simulated loops 
are shown in Figures 4.11 to 4.14. 


Table 4.1: Values of the constants in the constitutive equation 


(Constant name 

Symbol 

Unit 

Value 

mass 

m 

Kg 

2.2 

hcsteretic (lamping Co-efficient 

h 

Kg/.s^ 

2.11 X 10-'’ 

Conlornb friction Co-efficient 

fcol 

N 

1.5 


Cn 

N 

3.5 


C 2 

m~^s 

20.0 


Ck 

Tns~^ 

0.01 


Ca 

7n~^s 

90.0 


a 

Nim 

1476917.4 


b 

iV/m'-’ 

2839317.7 


9 


0.1 
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Figure 4.11: Experimental and simulated hysteresis loops ( 60 Hr, 0.1 mm deformation 
amplitude) 













Chapter 5 


Conclusions 


The dynamic characteristics of a rubber spring under harmonic loading have 
been studied. The shape of the hysteresis loops and the variation of the damping 
energy with deformation of the rubber spring leads to the following conclusions. 

[i] Both linear hysteresis damping and non-linear Coulomb damping, need to be 

combined for representing the damping force. The hysteresis dain|)ing was 
adopted in the model, since the energy dissipated per cycle failed to show 
any cliaracteristic variation with the frequency. Of course, tlie variation of 
damping energy with frequency showed considerable amount of scatter. 

[ii] The backbone curves of the hysteresis loops, while remaiiiing straight, changed 
orientation at different amplitudes. This behaviour has also been reported 
in earlier works. Conse(iuently, this amplitudes dependence of the dynamic 

stiffness needs to be modelled properly. 

[iii] A constitutive equation, iinder harmonic loading, has been proposed and 
seen to fit the result.s at room temperature within the. limitod range of fre- 
quency. 
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5.1 Scope of future work 


'I'lio proposod cxporiniental method should he rigorously tested in a wider range 
of frequency and a large number of samples which could not be done clue time 
constrant. The effects of temperature on the dynamic characteristics of rubber, 
known to be important, should be studied. Various configurations of standard 
metal-bonded rubber springs can be tested and proper models prepared for each 


configurations. 
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Appendix 


1. I. M. V. electrodynamic shaker 


Max. force output 

Max. acceleration 

Max. velocity 

Max. displacement 

Weight of the moving element 

Frequency range 

Max. loading weihgt 

Rating 

Power consumption 

2. Accelerometers 

Make 

Reference sensitivity 
Charge sensitivity 
Capacitance ( including cable 
Frequency range 


: 80 Kgf ( wihout load) 

: 47g 

: 41 cm/sec 

: 25 mm pp 

: 1.7 kg 

: 5-5000 Hz 
: 70 Kg 

: 24 Hours continuous duty 

( ambient temperature 0-40° C ) 
: 3.5 KVA 


B & K, Type 4370 

at 159.2 Hz, 100 ms“^ and 23° C 

10.20 pc/ms“^ 

1188 pF 
upto 25 kHz 
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Appendix 


*±U 


3. Oscilloscope 

Make 

Input impedence 

Max. allowable input voltage 

Cliannel section 

Time difference among two channels 
Sweep speed ( time base ) 

X-Y i)hase shift 


4. Charge amplifier 
Make 

Charge input 

Transducer sensitivity conditioning 
Amplifier sensitivity 

Calibrated output ratings 


Maximum output 
output impendence 


Kikusui COR 5501U 
1 mega-ohm ±2%,21pF±2pF 
400 V peak ( DC + AC peak ) 

CHI, ADD ( CHI ± CH2 ) 

CH2, CH2 INV, DUAL (CHOP, ALT) 
± 500 ps or less 
Real 20 ns/DIV - 0.5 s/DIV 
Storage 20 ns/DIV - 5 s/DIV 
DC-100 kHz within 3 degrees 


B & K 2635 

10® pc max. input rating 
3 digit dial in of transducer 
sensitivity from 0.1 to 0.99 pc/m.s~^ 
0.1 mV to 10 V/pC 
corresponding to -40 to 4-80 dB 
with transducer capacitance of 1 nF 
selectable in 10 dB steps 
Acceleration — o.l mV to 1 V/ms“^ 
Velocity — 10 mV to 100 
Displacement — o.l mV to 10 V/mm 
8 V ( 8 mA ) peak 
< 1 ohm • 










